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Accountancy and Power Flow Techniques for Aircraft Noise 
Transmission". The main emphasis of this report is on the use of 
the mobility power flow approach to deal with the fluid structure 
interaction problem. In conformity with previous reports, a case 
study of an L-shaped plate will be considered to demonstrate the 
approach. This report is for the period between June 1989 to 
December 1989. During this period one Master's thesis has been 
finished and will be defended the first week of January 1990. 
The topic of this finished thesis is the derivation of the 
mobility power flow approach for the case of acoustic excitation 
including scatter. A second Master's thesis which addresses the 
experimental side of the acoustic excitation problem is in the 
final stages and will be finished the first quarter of 1990. 
This is the seventh in this series of progress reports under this 
research grant. 
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ABSTRACT 


An analytical investigation based on the Mobility Power Flow 
method is presented for the determination of the vibrational 
response and power flow for two coupled flat plate structures in 
an L-shaped configuration, subjected to acoustical excitation. 
The principle of the mobility power flow method consists of 
dividing the global structure into a series of subsystems coupled 
together using mobility functions. Each separate subsystem is 
analyzed independently to determine the structural mobility 
functions for the junction and excitation locations. The 
mobility functions, together with the characteristics of the 
junction between the subsystems, are then used to determine the 
response of the global structure and the power flow. In the 
coupled plate structure considered here, mobility power flow 
expressions are derived for excitation by an incident acoustic 
plane wave. In this case the forces (acoustic pressures) acting 
on the structure are dependent on the response of the structure 
because of the scattered pressure component. The interaction 
between the structure and the fluid leads to the derivation of a 
corrected mode shape for the plates' normal surface velocity and 
also for the structure mobility functions. The determination of 
the scattered pressure components in the expressions for the 
power flow represents an additional component in the power flow 
balance for the source plate and the receiver plate. This 
component represents the radiated acoustical power from the plate 
structure . 



INTRODUCTION 


Mobility Power Flow (MPF) methods [1-4] have been shown 
previously to provide an effective structure analysis tool to 
deal with the response of one and two dimensional structures with 
direct mechanical excitation. Comparison of the MPF results with 
experimental results on a one dimensional structure [4], and with 
numerical (Finite Element Analysis (FEA)) [3] and statistical 
(Statistical Energy Analysis (SEA)) [4] results showed good 
agreement. Using the MPF method, the exchange of vibrational 
power between the substructures can be obtained for different 
structural wave components [5]. 

One limitation of past work [1-5] is that only excitation by 
direct point or distributed forces has been considered and that 
interaction with the surrounding medium has been neglected. In 
practical situations, distributed loads from acoustic excitation 
often apply on structures, where the influence of the surrounding 
medium may not be negligible. Therefore, the main objective of 
this study is to extend the MPF method to distributed excitation 
when the excitation is influenced by the structural response 
(acoustic excitation) . To be consistent with previous work the 
structure that will be considered in the analysis is an L-shaped 
plate, with an acoustic wave incident on one side of one of the 
flat plate substructures. 

In using the MPF, the structure is modeled by a series of 
coupled substructures with each substructure analyzed 
independently. The coupling between the substructures is 
defined through the boundary conditions at the junctions, taking 
into account the forces and moments that the substructures exert 
on each other. Expressions for the input power to a substructure 
and for the transferred power between the substructures, are 
obtained in terms of the input and junction velocity 
contributions and the forces and moments that are applied on the 
substructures . 

With direct force excitation of the structure, the input and 
transferred power can be written in terms of the input and 
transfer structural mobility functions of the substructures 
[2]. These only depend on the geometry of the specific 
substructure and on the frequency of excitation. The solution 
for the response of the global structure can be presented in a 
matrix form, which allows the method to deal with a large number 
of connected substructures. One advantage of dividing the global 
structure into subsystems is that, if the dependence of the 
response of the global structure on one of the subsystems is 
required, it is not necessary to repeat the whole of the analysis 
for the global structure, but only for that part that deals with 
the modified subsystem. Furthermore, the analysis is efficient 
to implement if the subsystem elements are identical, in which 
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case the independent response of only a typical subsystem is 
required. 

Having defined the approach for a distributed mechanical 
excitation in a past report [6], the analysis is extended here to 
deal with the acoustic excitation of the structure. A harmonic 
plane acoustic wave is assumed to be incident on one side of one 
of the flat plate substructures. With acoustic excitation, the 
effects of fluid loading may not be negligible. The scattered 
pressure will interfere with the incident pressure field and the 
mode shape of the structure will be influenced by the presence of 
the fluid. The fluid that will be considered here is air and 
therefore approximations applicable for light fluid loading can 
be applied. 


ACOUSTIC EXCITATION 

Consider a structure in the plane (x,y,0), with an acoustic 
fluid occupying the half-space z > 0. A plane wave incident on 
the surface of this structure can be described by: 


Pi (x,y , z) = P Q e(~ikxcos8sin<t>) e (-jkysinflsin^) e ( kzcos< £) 1. 


The direction of the plane wave incident on the surface of 
the structure subtends an angle <f> to the normal of the plate and 
an angle 8 to the x-axis direction in the plane of the structure 
(figure 1) . This acoustic wave, incident on the surface of the 
structure represents the excitation and the structure is set in 
motion. The motion of the structure creates an acoustic field 
which interact with the incident field. Thus, the total 
pressure acting on the structure surface is modified due to the 
presence of the scattered acoustic pressure component. The total 
pressure acting on the surface of the structure can be 
represented by [7]: 


p (x, y , z=0) = 2pi(x,y,z=0) + p s (x,y,z=0) 


where Pi(x,y,z=0) represents the incident pressure component and 
Ps( x »y/ z= 0) represents the scattered pressure component. The 
factor of 2 in front of the incident pressure component is 
introduced to take into account the reflected pressure component 
which is equal and opposite in the z direction to the incident 
pressure component (twice Pi(x,y,z=0) represents the blocked 
pressure component) . 

From equation (2) it can be observed that the scattered 
pressure component can be dealt with separate from the incident 
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pressure. If the structure has responses which are associated 
with different forces and moments acting on the structure, each 
of these responses creates a scattered pressure component. The 
total scattered pressure, thus has contributions from these 
different components and these can be dealt with separately. 

In the case of the L-shaped plate, the scattered pressure 
from the source plate, the plate receiving the acoustic wave, has 
two components, one which can be associated with the acoustic 
excitation and another component which can be associated with the 
edge moment representing the influence of the attached receiver 
plate. Each of these components can be dealt with separately, 
provided the changes in the response of the structure due to the 
fluid loading are taken into account. 

If the structure response mode shape is apriori known, the 
influence of the fluid loading on the structure can be directly 
determined from a solution of the coupled equations of motion 
[7]. However, this approach is not applicable when the mode 
function is not known. An alternative approach has been 
suggested by Leppington [8], where a solution for the response 
and scattered pressure are obtained based on an approximate 
solution for the case where the scattered pressure is not taken 
into account. The complete solution is given in terms of a 
correction factor, which is introduced to take into account the 
scattered pressure component and still satisfy the boundary 
conditions for the structure and the acoustic medium. 

With the Leppington [8] approach, the solutions for the 
structure response and scattered pressure are obtained by first 
determining the influence of the fluid loading on the response of 
the structure. This influence is described by a correction 
factor applied to the response of the structure obtained when the 
fluid loading is neglected. This correcting factor is 
independent of the magnitude of the applied loading and is mainly 
controlled by the boundary conditions for the structure and the 
type of loading. From knowledge of the response, a scattered 
pressure is computed. This scattered pressure is used to 
determine the pressure loading on the surface of the structure 
and also the vibrational power input to the structure. 

The above approach would directly fit into the MPF method, 
since with the MPF the global structure response is analyzed in 
terms of the separate responses of the subsystems representing 
the global structure. The scattered pressure components of each 
subsystem is thus separately determined and then summed to 
compute the power input to the global structure. That is, in 
this case of acoustic excitation of the L-shaped plate, the 
approach that will be followed would be to separately determine 
the response, including the correction factor, for each of the 
two plate subsystems, one representing a simply supported plate 
with an incident acoustic wave and the second representing a 
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simply supported plate with an edge moment. Then after having 
determined these corrected mobility functions, the total response 
of the structure and the scattered pressure components are 
computed. The power flow components are then determined from 
these results . The solutions for the two subsystems of the L- 
shaped plate follow in the next sections. 


SIMPLY SUPPORTED FLATE SUBJECTED TO AN INCIDENT SOUND FIELD 

The equation of motion for the transverse displacement of a 
plate is given by: 


D p VW(x,y) - pwW(x,y) = p(x,y, z=0) 3. 


where p(x,y,z=0) represents the pressure loading the plate. This 
term consists of the incident pressure component plus the 
scattered pressure component (equation 2) . 

Assuming a sinusoidal mode shape in the x direction [2], the 
direction parallel to the junction, . the displacement W(x,y) can 
be written in the form. 


W(x,y) 


W (y) sin 
m J 


r 'v 
mTTX 


m=l 


Substituting this equation into equation (3) , 


r i 

4 

f y 
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m ir 

W(y) - 2 

mn 





. a . 


. a , 



9 W(y) d W(y) p(y) 

+ k W(y) 


3y 


3y 
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where the subscript m is dropped for simplicity, kp is the plate 
wavenumber and p(y) is defined by 


p(y) - 
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r a 

p(x,y)sin 


mTrx 


a J 


* o 


dx 
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Applying equation (6) to the incident pressure, given by equation 

( 1 ), 


p t (y) - 


?p 

o 


a 


t (-jkysin0sin<£) 
L 1 e 


7. 


where 


nwr 


^ e ( -jkacos^sin^) 


(kcostf sin< 




- ( imr/a 
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The scattered pressure component for mode m is a function of 
the normal surface velocity and is therefore dependent on the 
solution for the response of the plate. That is, a coupled 
solution for the plate and the scattered pressure is required. If 
the scattered pressure component is neglected, the right hand- 
side of equation (5) has only one component, given by equation 
(7) multiplied by a factor of 2. In this case the solution for 
the displacement of the plate can be written in the form: 


W(y) = W g (y) + W h (y) 


9. 


where Wg(y) is the particular solution of equation (5) , 


w g(y > 


4P I, 

o 1 


( - jkysinflsin^) 


(mjr/a)^ + (ksintfsin^)^] ^ 


10 . 


and Wft(y) is the solution of the homogeneous equation of motion, 


W h (y) = Acosh(k^y) + Bsinh(k^y) +Ccos(k 2 y) + Dsin(k 2 y) 11. 


where k^and k 2 are defined by 



The coefficients A,B,C and D in equation (11) are selected to 
satisfy the boundary conditions in the y direction. In the case 
of the simply supported plate these boundary conditions are: 


W(y-O) - 0 


a 2 W(y=0) 



- 0 



13. 


Solving for A,B,C and D using these boundary conditions will 
yield 
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A - 


W''(°) + W (0)k* 

R R 2 


2k 


-k^W (b) - W' ' (b) + 
Jj g 


2l 

W' '(0) - W (O)k^lcosh(k-b) 
_g g 2| 1 


W' ' (0) 
g 


2k sinh(k 1 b) 
P 1 


W (0)k- 

& L 


c - 


2 k 


■kW (b) + W' ' (b) - 
1 5 I 


W' ' (0) - W (0)k, 
R R 1 


cos(k 2 b) 


2k sin(k„b) 
P 2 


14. 


where the double prime indicates a second derivative with respect 
to y. 

To compensate for neglecting the scattered pressure 
component, a corrective complex coefficient Kj is introduced [8] 


Wi(y) = W g (y) + KjWh (y) 


15. 


The main influence of the scattered pressure component will be on 
the contribution from the solution to the homogenous equation of 
motion [8], This is the reason why the correction factor K is 
only applied to the term Wjj(y). This complex correction 
coefficient K will be close to unity away from the resonant 
frequencies and will deviate from unity at the resonant 
frequencies where the influence of the fluid on the velocity mode 
shape is most significant [8]. If W^(y) represents the exact 
solution for the response of the simply supported plate to the 
acoustic excitation, then the approximate result obtained when 
neglecting the scattered pressure component introduces an error 
w e(Y) gi ve n by 


W e (y) = W(y) - Wi(y) 


16. 
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By substituting the above equation into the equation of 
motion (equation 5) , the error term W e (y) must satisfy the 
differential equation, 


^4 


nwr 

W e (y) - 2 

r ^ 

mn 
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w a - 


2 2 

cTW e (y) 



+ 


a 4 w e (y) 



k 4 V (y) 
p e 


p s ( y) 


D 

P 

17. 


In solving equation (17) , a set of boundary conditions that 
apply to W e (y) can be derived from the boundary conditions for 
W(y) • 


W (y=0 or b) 
e 


W 1 (y-0 or b) - -(W (0 or b) + KjW^O orb)) - 0 


dVT(y=0 or b) 


ay 


aw^(y=0 or b) 


ay 


aw (y=0 or b) 

e 


K, 


aW h (y-0 or b) 


ay 


ay 


18, 


A solution for W e (y) is derived in terms of the Green function 
G(£,y) for the structure [8] 


W e (y) 


1 


D 

P 


b 

G((,y)p (Od( + G ' (b ,y)W' ' (b) - G ' (0 ,y)W' ' (0) 

S 6 c 

o 

+ G" ' (b,y)W" (b) - G' " (0,y)W" (0) 
- 2 (nwr/a) |g ' (b , y)W e (b) - G' (0,y)W e (0) j 

19. 


where the primes for the function G represent derivatives with 
respect to £ . The function G (£,y) is defined by: 
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From the expression for W e (y), the coefficient K is selected 
such that W e (y) = 0 for all values of y. Substituting equation 
(20) into equation (19), and equating W e (y) to zero, an 
expression is obtained that contains p s ((), which is still an 
unknown. However, p s (£) can be expressed as a function of the 
normal displacement of the plate, by using the momentum equation. 
Introducing this substitution and after some manipulations which 
are presented in the appendix, the following result is obtained 
for p s (£) 


1 


D 

P 


r b 

P s (0sin(k 2 0 
J o 


= 


2n 2 a 



W 1 (^)I(/3)S(^)d^ 



where is the fluid loading coefficient defined by: 


a l = (wp o)/ D p 


22 . 


p Q is the fluid density, and p is the spatial Fourier transform 
variable in the Y direction. I(/9) and S(/3) are defined in the 
appendix. From equations (15), (18), (19) and (21) an expression 
for the correcting coefficient Kj is derived. 
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J± 


f ^ 

m7r 


2* 2 a 


W (P)I(p)S(0)d0 + RES 


K, 


2?r 2 a 


ITITT 


1 2 


V J 


W h (^)I(^)S(^)d^ + RES 


23 . 


where 


RES 


k 0 cos(k 9 b)W' '(b) - kH-'(O) - k 2 k cos(k b)W (b) + k 2 k W (0) 
2 2s 2 e 12 2 e 1 2 e 


24 . 


Having solved for the response of the simple supported plate 
a modal mobility function can be defined by the ratio of the 
modal velocity response to the input incident pressure. That is, 
using the same notation as in previous reports, with the 
subscript 1 indicating the input excitation location, the input 
modal mobility for mode m is given by. 


(y) - 

m 


jwW^y) 


P 

o 


jw4I 1 e 


(-jk a y) 


aD 


+ jwK 


A'cosh(k^y) + B'sinh(k^y) + 
C'cos(k 2 y) + D'sin(k 2 y) j 


25 . 


where 
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k - ksin0sin<4 
a 


26 . 


2 I, 


aD k 2 
P P 


2 2 

1c - V 
2 a 


27 . 


21 , 


B' 


2 

aD k 
P P 


- 21 , 


C' 


21 , 


D' 


2 ^2^ Q -jk a b 
a 2 ‘ 


(k - k„)e - (k^ + k^)cosh(k 1 b) 

" a 2 l 


3 sinh(k^b) 


2 

aD k 
P P 


2 2 

kf + k 

1 a 
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29. 
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aD k 
P P 


(k 2 + k 2 )e + (k 2 + k 2 ) cos(k 0 b) 

a 1 a 1 2 


sin(k 2 b) 


30 . 
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ITITT 
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+ k 
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a 
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Similarly, a modal transfer mobility function, defined by the 
ratio of the edge rotational velocity per unit incident pressure, 
can be derived from the above solution for the response. If the 
subscript 2 represents the location of the junction, the modal 
transfer mobility for mode m is given by. 
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aw 1 (y-b) 


J w 


M 


21 


m 


dy 


w4I, 


aD 


k e 
a 


-jk a b 


jw Kj A'k^sinh(k^b) 


B 'k^cosh(k^b) - C'k2sin(k2b) + D'k2Cos(k2b) 


32 . 


One can observe that for these two mobility functions the 
input mobility is a function of both the mode number and the 
variable y while the transfer mobility is a function of only the 
mode number. The reason for the y dependency for the input 
mobility is because the excitation is distributed over the 
surface of the plate structure and thus this input mobility 
represents the response per unit incident pressure anywhere on 
the surface of the plate structure. 


SIMPLY SUPPORTED PLATE WITH AN EDGE MOMENT 

An approach similar to the one used for the subsystem 
discussed in the previous section will be used for this second 
subsystem. In this case there is no incident acoustic waves, 
that is Pi(x,y,z) = 0. The forces acting on this subsystem are 
the edge moment and the scattered pressure. Therefore the same 
eguation of motion as equation (5) applies with the pressure term 
on the right hand side only representing the scattered pressure. 
Because of the presence of the edge moment, the boundary 
conditions for the edge y=b of this plate subsystem are different 
from those given in equation (13). The boundary conditions for 
the edge y=b in this case are given by. 


a 2 W(y=b) T 

W(y-b) - 0 ; 

2 

dy D 

P 


33 . 


where T m represents the mode m component of the edge moment [2]. 
the response in this case is obtained by solving the equation of 
motion as given by equation (5), with p(x,y,z=0) replaced by the 
scattered pressure component. In obtaining a solution to this 
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equation of motion, a first approximation is obtained by 
neglecting the scattered pressure, in which case it becomes a 
homogeneous equation of motion. The solution that satisfies the 
boundary conditions is given by, 


T 

W h (y) = W(y) m 

2J pD w 
P 


sin(k 2 y) sinhCk^y) 


sin(k 2 b) sinh(k^b) 


34. 


To account for the presence of the scattered pressure 
component, a corrective factor Kjj is introduced 


w L (y) 


kt 

II m 


2J pT) u> 
P 


s..n(k 2 y) sinh(k^y) 


sin(k 2 b) sinh(k^b) 


35. 


Wi(y) represents the exact solution of the equation of motion 
including the scattered pressure. The error introduced by the 
approximation is again given by an equation of the form of 
equation (16) where in this case as well W e (y) must satisfy 
equation (17) . However, the boundary conditions that apply to 
the term W e (y) in this case are different from those given in 
equation (18) . The boundary conditions are modified to include 
the influence of the edge moment. 


W g (y=0 or b) - - W^y-0 or b) - - K W (0 or b) - 0 


dW^(y-O) 


dW^(y-O) 


dy 


dy 


3W (y-b) 
e 


dW^(y-b) 


dy 


dy 


aw. <y=0) 
0 

dy 2 


m 


m 


a - K n ) 


36. 


The solution for W e (y) is again derived in terms of the 
Green function (similar equations as (19) and (20)). The 
application of the boundary conditions into equation (19), and 
setting W e (y) = 0 for all values of y leads to an equation 

involving p s (y), 
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D 

P 


f b 

P s (Osin(k 2 C)d| - k 2 cos(k 2 b)(l - 
o 


0 


37. 


Substituting for the integral in equation (37) by 
equation(21) , an expression for the correction coefficient Kjj is 
obtained as follows, 


K 


II 


k 2 cos(k 2 b) 


k 2 cos(k 2 b) + 


j 


2tt 2 a 



W h (^)I(^)S(^)dy3 
- 00 


38. 


Having derived a solution for the response of the plate 
subsystem sub j ected to an edge moment and including the influence 
of the fluid loading on the plate, an input mobility function can 
obtained for the edge of the plate. Defining an edge input 
mobility as the rotational velocity response per unit applied 
edge moment, 


M 


2 


m 


jw dW(y=-b) 

T dy 
m 





tan(k 2 b) tanh(k^b) 


39. 


If the two plates of the L-shaped plate structure are identical 
then. 



m m 


40. 


The subscript 3 represents the connected edge of the receiver 
plate [2 ] . 
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A transfer mobility function for the plate surface velocity 
response per unit applied edge moment can also be defined from 
the above analysis, 


M 


12 


m 



sin(k 2 y) sinhCk^y) 


sinCkjb) sinh(R^b) 


41. 


In this case as well the input mobility functions are only a 
function of the mode number m, However, the transfer mobility is 
a function of both the variable y and the mode number m, since it 
represents a transfer to any point on the surface of the plate. 


INPUT AND TRANSFER POWER EQUATIONS 

Having derived the mobility functions, the derivation of the 
power flow equations follows in the same way as was done for the 
mechanical excitation [2]. There are however some differences. 
One of the differences is that some of the mobility functions, 
apart from the modal dependency, also include a dependency on 
the spatial variable y. Since the excitation is a distributed 
load the input power must be obtained both as an integral over 
the spatial variable and as a summation over all the modes. The 
modal summation is an alternative way of performing a spatial 
integration when the response can be decomposed into a set of 
modes. 

For the transferred power, this is given by an integral along 
the entire length of the junction or alternatively since this is 
the direction for which a modal decomposition has been assumed, 
as a summation over all the modes. Since the transferred power 
is dependent on the edge moment which is controlled by the 
incident acoustic excitation, to evaluate the edge moment and 
hence the transferred power an integral still has to be performed 
for the y direction, the direction perpendicular to the junction. 

Power Input 

The total input power is given by the product of the total 
pressure acting on the source plate surface and the plate 
velocity response integrated over the y direction and summed for 
all modes m. 
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Real 
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2p i (y) + p s ^(y) + p s ^(y) ] v *(y) d y 

^ o 

42 . 


where the two p s terms represent the scattered pressure 
components associated one with the response of the source plate 
due to the incident sound wave and one due to the application of 
the edge moment. Vj* is the surface velocity of the source plate 
for mode m and is given by, 


Vy> - p o M i <y) + Vu w 43 • 

m m 

Eliminating Tm by solving for continuity of motion at the 
junction edge, 
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m 



T M„ 
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P M 
o 21 


m 


2 M„ 
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44 . 


substituting into equation (43) , 


v 1 (y) = P 


o 


M 1 (y) + 
m 



m 



M 12 (y) 

m 


45 . 


Substituting into the power input expressions the values for 
Vi (y) and for the scattered pressure components the following 
result is obtained for the input power; 
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46. 

In the last integral of the above equation the variable y has 
been changed to the variable p which represents the fourier 
transform variable with respect to y. Also, this last integral 
represents the power flow out of the source plate, as radiated 
acoustical power. 


Power Transfer 


The transferred power between the two plates is given by the 
summation of all the mode contributions for the product between 
the edge moment and the edge rotational velocity. This can also 
be expressed in terms of the modal mobility functions derived 
above , 
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trans 


a 

4 


00 




O 




2 

Real 
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RESULTS 

Results for the power input, power transfer and radiated 
acoustical power are obtained using the above power flow 
expressions. The characteristics of the plate structure 
considered in the analysis are the same as those for previous 
analysis in this series. These characteristics are shown in 
figure (2. 
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Acoustic waves are considered incident on the source plate 
such that both 9 and <f> are equal to 45°. The results for the 
input and transferred power for this angle of incidence are shown 
in figure (3). Figure (4) shows the acoustic power radiated by 
the source plate obtained from the scattered pressure component. 
As can be observed from this figure the contribution to the 
scattered pressure from motion of the source plate due to the 
presence of the edge moment is not significant. 

To determine the influence of the fluid loading on the power 
flow, the input and transferred power are computed when the the 
fluid loading effects (both the scattered pressure and the 
coefficients K) are neglected. Figures (5) and (6) show a 
comparison between the two sets of results. As can be observed 
from these figures the influence of the fluid loading is not very 
significant, although some differences are observed mainly near 
the resonant frequencies and in the troughs between the resonant 
frequencies. The main influence of the fluid loading is to 
increase the damping of the structure due to the acoustic 
radiation. 

Also considered are other angles of incidence. Figures (7) , 
(8) and (9) show the power input and power transfer for different 
values of the incidence angles 9 and <f>. The angles of incidence 
considered have values of 9 and <f> given by O'and 0°, 15°and 30°, 
and 75°and 60° respectively for figures (7), (8) and (9). The 
number of modes in figure (7) is reduced as compared to the other 
results mainly because of the symmetry of the excitation. The 
even modes are not excited with normal incidence. Apart from the 
number of modes that are excited the general shape of the power 
flow curves are also influenced by the angles of incidence. 


CONCLUSION 

The extension of the Power Flow Method to excitation 
conditions other than mechanical excitation has been demonstrated 
in this report. The excitation considered here is an incident 
acoustic wave. In this case of acoustic excitation, the response 
of the structure influences the incident acoustic field and the 
problem of the structure response becomes a fluid-structure 
interaction problem. 

Since most of the work found in the literature on fluid- 
structure interaction deals with simply supported plate 
structures, because of the requirement of an apriori knowledge of 
the vibration mode shape, which is not valid in the case of 
connected plate structures, an approximate solution based on the 
work by Leppington [8] has been used here. Although fluid 
loading is considered, it is found that light fluid loading does 
not significantly modify the mode shape of the vibrating plate 
structure but that the scattered pressure can be significant. 
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If the results obtained here for the power input and 
transferred for the L-shaped plate acoustic excitation are 
compared to results for mechanical excitation [6], the following 
observation can be made. With uniformly distributed mechanical 
excitation, the power flow is similar to chat for excitation from 
normal incidence acoustic waves. For oblique incidence waves 
additional modes of vibration are excited same as in the case of 
point excitation. Compared with the power flow results for point 
excitation, the power flow is higher in the case of the 
mechanical excitation. For these comparisons the total load on 
the source plate is kept constant. 

As a final conclusion, an important result of this report is 
that the effects of fluid loading on a connected plate structure 
can be integrated into the general Power Flow Method. This 
enhances the usefulness of power flow methods in determining the 
response of plate-like structural components and their 
interaction with the surrounding medium. 
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Appendix A 

RELATIONSHIP BETWEEN SCATTERED PRESSURE 
AND NORMAL SURFACE VELOCITY 


In deriving the relationship between pressure and normal 
surface velocity the following conditions are assumed (a) 
the plate normal displacement has a time dependency of the form 

W(x,y,t) = W(x, y) exp ( jut) (A.l) 

(b) the plate is finite and lies in the plane (x,y,0) with 
the fluid occupying the half-space z a 0; (c) the 

propagation of a plane acoustic wave in three dimensional 
space is given by: 

P (x; y , z , t) = pexp(-jkx -jk y -jk z) exp( jut) (A. 2) 

x y z 

The Fourier transform of W(y) is defined by: 

+ 00 

W(/3) * Jw(y)exp(j0y)dy (A. 3) 

- CO 


and the inverse transform by: 

+ 00 

W( y) - ^ Jw(0)exp(-j£y)d0 


(A. 4) 


i 


- 00 


From the momentum equation in the z direction: 



3 V 

at 


o 


(A. 5) 


where p is the acoustic pressure, p the fluid mean density 
and V the particle velocity in the z direction. From 
equation (A. 5) 

-jk z p s (x,y,z)exp(ju)t) = -j(jp Q V(x,y, z) exp( jwt) (A. 6) 

where p s is the scattered pressure. 

At the interface between the plate and the acoustic 
medium (z = 0), the surface normal velocity is equal to the 
acoustic particle velocity. Therefore at this interface 

up 

P s (X/Y) = ■ V(x,y) (A. 7) 

Z 

where k^ must satisfy the condition: 

k 2 + k 2 + k 2 = (£) 2 = k 2 (A. 8) 

where k is the acoustic wavenumber, <j the pulsation 
frequency and c the acoustic wave speed. 

Let a and /3 be the plate wavenumbers in the X and Y 
directions respectively. Since the acoustic wave is 
generated by the plate motion, the X and Y variations of 
the acoustic field must follow those of the plate and 


I 



therefore: 


k = a and k - j3 (A. 9) 

That is k 2 = k 2 - a 2 - /9 2 = k 2 - k 2 (A. 10) 

z p 

where k p is the plate bending wavenumber 

Spatial Fourier transforming equation (A. 7) in the X and 
Y directions 

“P 0 

P s ( a /&) ~ % (A. 11) 


For a finite rectangular plate, simply supported at 
x = 0 and x = a, the response of the plate can be described 
by: 


00 

V(x,y) = l V(y)sin(^) 
p = i a 


(A. 12) 


Fourier transforming the expression for V(x,y) 


V(a,<3) - f V(ff)(E2) ~ 1 

pti a a 2 - gZ .) 2 


For p (x,y), this is first written in the form 

s 


P s (x,y) = I P s (y)sin(Sip) 

m = 1 


(A. 13) 


(A. 14) 


l 



where 


P S ( Y) - | } p a (x,y)sin(5^)dx 


and then Fourier transforming in the Y direction 


P.(*> 


| P s (x,f3)sin(255)dx 


(A. 15) 


For the X direction, also taking a Fourier Transform 


a + co 

(£) = I j [ Ijt I P g (a,/3)exp(-jax) 

o L _oo 


da 


si n(25* )dx 

a 


(A. 16) 


-^{ P 5(a , g)( lE, I^E^|I -l . da 

-oo a ( a ' 

(A. 17) 

Combining equations (A. 11), (A. 13) and (A. 17) yields: 


PAP) = 


UQ 00 

, ItlTT v y 

na ^a ' ^ 

p=i 


(|E)V(0)* 


I 


(-l) m exp(jaa) - 
a 2 - (25) 2 


1 (-1) p exp(-jaa) 

2 - (E5) 2 
v a ' 


-11 


da 


a 


- 00 


The terms of the summation for which p * ra represent 
intramodal coupling which has been shown to be negligible 
for light fluid loading provided the modal density of the 
structure is low. 

Therefore equation (A. 18) can be simplified to 


up 

p.<*> - *ir<ir> 2l <*> v (e> 


(A. 19) 


where 


+ 00 

I(/3) = 2 | 
- 00 



m 


(-1) cosaa ] 


- 0 3 



(A. 20) 


From (A. 20) equation ( 21 ) can be derived since the 

definition of the inverse Fourier transform is: 

+ 00 

P.(*> = hi Jp s 0) e xp(-j/3?)d|3 (A. 21) 

- 00 

and substituting for p # ($) 
b 

|p s (C)sin(k 2 ^)dC = 

0 

+® b 

h iFirO 2 |lO) v 03) [ Jsin(k 2 Oexp(-j/3C)dc]d(S 
— 00 0 


(A. 22) 



or, since V(/3) = ja>W(/9) , where W is the displacement 


1 

D 


b 

|p s (C)sin(k 2 e)de = 

0 


2 . 

j__ (jnTT) 2 (___£) Ji(/3)W(/3)S(i3)d/3 


2rr a 


(A. 23) 


where 


SO) = Jsin(k 2 ^)exp(-j/3C)d? 




- k 2 + exp(-j/9b) [j/3sin(k b) + k cos(k b)] 


(A. 24) 


It should be noted that S(/3) is not the Fourier 
transform of sin(k 2 £) . From this expression the scattered 
pressure component can be evaluated. 




Materiel: Aluminium 
Density: 2710 Kg/m » 

Elastic Modulus: 7.2 x10 1# N/m 1 
Thickness: 0.00635m 
Dimensions: a«1.0m, b«=0.5m 
Loss F actor. 0.01 

Figure 2. Plate structure showing plate characteristics, 
and coordinate directions. 
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Figure 4. Radlittd scout foal power per unit eoouetio preeeure by the 

eouree piste. I Totsl redieted power i s contribution to 

rsdlsted power frow edge wo went . 






Figure 9. Power Input per unit scouetle pressure with end without the 
Influence of fluid loedlng, » With Fluid loedlng effeotei — -I 

without fluid loedlng effects. 



Figure 6, Power trenefer per unit eoouetlo pressure with end without 

the influence of fluid loedlng* i With Fluid loedfngi — i without 

fluid loedlng* 
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Figure 7 , Power Input and traefer per unit aoouetlo preeeure th th 
Influence of fluid loading for angle of fnoldenoe 0 and 9 equal 0*and 
0 * reepeotlvely. i Power input i - — i power transfer. 
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Figure 8. Power Input and traefer per unit aoouetlo preeeure with the 
Influence of fluid loading for angle of inoldenoe 0 and 6 equal 15 * and 
30 # reepeotlvely. i Power Input 9 1 power tranefer. 
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FI gura 9, Poiar Input and traafar par unit aoouatfo praaaura alth tha 
influanoa of fluid loading for angla of tnoldanoa ft and 9 aqua! 79* and 
60* raapaot f va 1 y • t Pouar Input i — i powar tranafar. 





